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Abstract. Constraints on the dynamics, dissipation, and production of VHE 
neutrinos in relativistic jets are derived using opacity calculations and VHE 
7-ray observations. In particular, it is demonstrated how rapid variability of 
the 7-ray emission at very high energies (> 100 GeV) can be used to map the 
location of the 7-spheres, to derive lower limits on the Doppler factor of the 7-ray 
emission zone, and to constrain the photopion production opacity. The apparent 
discrepancy between jet Lorentz factors inferred from superluminal motions and 
source statistics in the TeV blazars and those derived from the 7-ray emission 
is discussed. The relation to the high-energy emission from the HSTl knot in 
M87 is briefly mentioned. Estimates of neutrino yields in upcoming neutrino 
telescopes are given for various sources. It is shown that for TeV blazars, the 
rapid variability of the TeV emission implies neutrino yields well below detection 
limit. 



1. Introduction 

There is now strong evidence that the high-energy emission observed in blazars, 
microquasars, and 7-ray bursts originates from collimated, relativistic outflows. 
The common view is that those outflows are powered by a magnetized accretion 
disk and a spinning black hole and collimated by magnetic fields and/or the 
medium surrounding the jet. However, there is as yet no universal agreement 
about the mechanisms responsible for the formation and acceleration the jet and 
the dissipation of its bulk energy. Even the composition of the jet is unknown in 
most sources. The rapid variability frequently observed suggests that the high- 
energy emission is likely to be produced close to the inner engine and, therefore, 
provides an important probe of the physics involved on those scales. As discussed 
below, some important constraints on the jet parameters can already be imposed 
using available data, and further progress is expected in the coming few years 
with the advance of observational techniques. A new generation of experiments 
just started operating or will become operative soon: Space-based (e.g. GLAST, 
INTEGRAL, AGILE), Cerenkov (e.g. HESS, MAGIC, VERITAS), and air- 
shower (e.g. MILAGRO) 7-ray detectors will probe with high sensitivity the 
energy range of 10 MeV to a few TeV (see, e.g., Catanese, 1999 for a review). The 
operating BAIKAL and AMANDA neutrino telescopes, and the cubic-km scale 
telescopes under construction (IceCube, ANTARES, NESTOR, NEMO; see, e.g. 
Halzen, 2005 for review), will open a new window onto the Universe. Besides 
providing an important probe of the innermost regions of compact astrophysical 
systems, these experiments can also be exploited to test new physics. Finally, 



1 



2 



Amir Levinson 



new ultra-high-energy cosmic-ray detectors, e.g. the HiRes and the hybrid Auger 
detectors will provide cosmic-ray data of unprecedented quality and quantity. 

2. Gamma Rays 

At sufficiently small radii the compactness parameter, = L^a^ / {rn(,c^r) ~ 
10^(L-y/LE(jd)(''/''^s)~^) largely exceeds unity for essentially all classes of compact 
high-energy sources. It is, therefore, expected that 7-rays will not be able to 
escape from the inner jet regions without creating pairs. Both the synchrotron 
photons produced inside the jet and the ambient radiation intercepted by the jet 
contribute an opacity to pair production. The threshold 7-ray energy for pair 
production via interaction with a target photon of energy hvg is 

e,,. - 2.5 X 10^1 (^)"' eV. (1) 

Consequently, for a target photon spectrum Fi, oc u~'^ with a > the 7-ray 
opacity at a given radius should increase with increasing 7-ray energy. Results 
of detailed calculations of the 7-ray opacity are exhibited in figure 1, where the 
^-spheric radius, defined as the radius r^(e^,) beyond which the pair production 
optical depth to infinity is unity, viz., r-y-y(r-y, e-y) = 1, is plotted against 7-ray 
energy e-y, for two target radiation fields: synchrotron radiation (dashed lines) 
and external radiation (solid lines). The spectra of the target radiation fields 
employed in those calculations are discussed in Levinson (2006). As seen the 
7-spheric radius indeed increases with increasing 7-ray energy. 

The 7-spheres can be mapped in principle by measuring temporal variations 
of the 7-ray flux in different energy bands during a flare. If the 7-ray emission is 
produced over many octaves of jet radius, where intense pair cascades at the ob- 
served energies are important (Blandford & Levinson 1995), then it is expected 
that a flare will propagate from low to high energies, or that the variations at 
higher 7-ray energy will be slower than at lower energies. With the limited sen- 
sitivity and energy band of the EGRET instrTimcnt it was practically impossible 
to resolve such effects. It is hoped that with the upcoming GLAST instrument 
this will be feasible 

Further constraints on the source parameters can be derived by measur- 
ing the low-energy flux (radio-to-IR) simultaneously with the variable 7-ray 
emission. Suppose that a variability timescale Ai has been measured at some 
observed 7-ray energy e^. This implies that the emission at the observed energy 
originated from a region of size (as measured in the Lab frame) Ar < ryar = 
T5c/S.t/ {1 + z), where T and 5 are the bulk Lorentz factor and the corresponding 
Doppler factor of the emitting matter, respectively, and z is the redshift of the 
source. Above the emission region the pair production optical depth at the ob- 
served energy must not exceed unity, viz., T^^'^(rem) e^/^) < 1, where rem is the 
radius at which the 7-ray emitting region is located. Assuming rem — Ar < rvar 
implies r^(e'y) < ryar- Suppose now that the synchrotron flux near the peak can 
be measured simultaneously. The latter condition on r-y(e-y) can be solved for 
the Doppler factor to yield (Levinson 2006) 

55 > 2 X lQ^^Atl^''^{Te)-^z^{e^/lTeWfl^Siy, (2) 
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Figure 1. Dimcnsionlcss 7-sphcric radius versus 7-ray energy, computed 
using the target photon spectra from Levinson (2006). The different curves 
correspond to a different normaUzation of the target radiation field intensity. 

where 5jy is the measured synchrotron flux density in Janskys, 9 is the opening 
angle of the jet and Ath is At measured in hours. Adopting ~ 1, we estimate 
(5 > 35 for the rapid flare observed in Mrk 421 and 5 > 190 for the few minuets 
variability reported for PKS 2155-304, with rem < 10^'' cm for the minimum 
condition in both sources. Such high values of S are in clear disagreement with 
the much lower values inferred from source statistics (e.g., Urry & Padovani 
1991) and superluminal motions on parsec scales (Jorstad et al. 2001). Various 
explanations, including jet deceleration (Gcorganopouloselta & Kazanas 2003), 
a structure consisting of interacting spine and sheath (Ghiscllini ct al. 2005), 
and opening angle effects (Gopal-Krishna ct al. 2004) have been proposed in 
order to resolve this discrepancy. Alternatively, it is conceivable that those flares 
result from impulsive ejection episodes of high Lorentz factor shells, like in the 
internal shock model of GRBs. 

It has been argued that such high values of the Doppler factor may not 
be required if the 7-ray production region is located far from the black hole, 
at radii rem >> Ar. However, the fraction of jet energy that can be tapped 
for production of 7 rays in a region of size Ar located at a radius rem is rj ^ 
{Ar/Ore.mf'- As a consequence, either the opening angle of the jet is very small, 
6 ~ Ar/r << 1 or the jet power is much larger than the luminosity of the 
TeV emission measured during the flaring states. Another possibility is that the 
TeV emission is produced by a converging shock in a reconfinement nozzle, as 
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proposed for the HSTl knot in M87 (Cheung et al. 2007). It should be noted that 
in M87 the X-ray and TeV luminosities are -LxeV ~ -^^x ~ 10^^ erg s~^ (Aharonian 
et al. 2006; Cheung et al. 2007), much smaller than the TeV luminosity LxeV ~ 
-|^q44-45 gj.g g-i oi^ggrved typically in the TeV blazars. Estimates of the jet power 
in M87 yield Lj ^ 10^^ erg s~^ (e.g., Stawarz et al. 2006; Bicknell Sz Begelman 
1996), implying a very small conversion fraction, L^eV / Lj ^ W~^. Even with 
such a small conversion efficiency an opening angle 9 < 10^^ rad is required if 
the TeV emission were to originate from the HSTl knot, unless reconfinement, 
as mentioned above, can give rise to convergence of the jet. This may occur if 
radiative cooling of the shocked jet layer is effective. Recollimation shocks may 
be an important dissipation channel also in other sources, e.g., GRBs (Bromberg 
& Levinson, 2007) 

3. Cosmic rays and Neutrinos 

A substantial fraction of the energy dissipated in the jet can, in principle, be 
tapped for the acceleration of protons to ultrahigh energies with a power law 
spectrum. The maximum proton energy is limited by energy losses at small radii 
and by escape at larger radii (Levinson 2006), and in any case is restricted to 
ep,max < eBOr, where 9 is the opening angle of the jet and B the magnetic field. 
In terms of the equipartition parameter the magnetic field can be expressed as 
B = AxW^{^BCj/mBH)^^'^{9r)-^ Gauss, where Cj is the jet power in Eddington 
units and niBH is the mass of the black hole in units of solar mass. Consequently, 
production of UHECRs of observed energy requires 

^BniBHCj > 10«(ep/1020-5eV)2, (3) 

which essentially leaves GRBs with uibh ~ 3, ~ 0.1, Cj ~ 10^^, and pow- 
erful blazars with tubh ~ 10^, — 0.1, Cj ~ 1 as the main candidates for 
astrophysical UHECRs source^ Whether protons can actually be accelerated 
by relativistic shocks to the highest energies observed is yet an open issue. 

High-energy neutrinos can be produced in astrophysical jets mainly through 
the decay of charged pions. The pions may be produced through collisions of 
protons with target photons (p'j), or via inelastic {pp) and (pn) collisions. The 
former process is usually dominant in jets (for some exceptions see, Levinson &: 
Eichler, 2003; Torres & Halzen, 2007). 

Observations of VHE 7-rays can constrain the photopion production opac- 
ity, particularly in situations where rapid variability of the VHE 7-ray flux is 
observed. Because both the protons and the 7 rays interact locally with the same 
target radiation field the ratio of photomeson and pair production opacities de- 
pend solely on the ratio of cross sections, cTpj/a^^ ~ 4 x 10^^, and the spectrum 
of the target radiation field. For a target photon spectrum ns{v) oc this 
ratio reads 

T^^{ey,r) ns{e^,th) (y-y-y \^3 x lO^e^ 




^ Condition Q does not apply to vacuum gaps in dormant AGNs (Boldt & Gosh 1999; Levinson 
2000) and magnetars where the ideal MHD condition is violated. 
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Detailed calculations of opacity ratios that employed more realistic target photon 
spectra are presented in Levinson (2006), and are plotted in fig. 2. As seen from 
the figure, at 7-ray energies above a few TeV the opacity ratio is smaller than 
unity even at the maximum proton energy. For the class of TeV blazars this 
implies neutrino yields well below detection limit (see Levinson 2006 for more 
details) . With GLAST it should be possible to constrain other sources and to use 
such constraints to identify the best candidates for the upcoming km^ detectors. 




Figure 2. Opacity ratio, Tp^^ (e'p) / t^^^" {e'^) , as a function of comoving pro- 
ton energy e^' at two radii: r/r^ = 10^ (dashed lines) and 10^ (solid lines). 
The left panel corresponds to a choice of parameters that represents a typical 
blazar like 3c279, and the right panel to a choice of parameters representing 
a typical microquasar like GRS 1915. The numbers that label the curves 
correspond to the comoving 7-ray energy in units of TeV. The dotted Unes 
designate the maximum proton energy derived from the confinement limit. 



The neutrino flux emitted from the jet depends on several factors: (i) 
the fraction of jet energy injected as a power law distribution of protons ^p, 
(ii) the spectrum of accelerated protons, (iii) the cooling times of pions and 
muons and (iv) the pion production opacity. It is customary to introduce 
the parameter fniep) < 1, representing the fraction of proton energy con- 
verted to pions. In case of photopion production it can be approximated by 
f-^i^'p) = min[l, KT^Tp^^e'p)], where K.,^ is the elasticity. For a proton distribution 
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np{ep) the average fraction of proton energy lost to pions can be defined as 

This average fraction depends on the spectrum of both protons and target pho- 
tons. For powerful blazars this fraction is estimated to be /tt ^ 0.2 assuming 
a proton energy distribution np{e) cx e~^, whereas for the TcV BL Lac sources 
/tt < lO^"^ (Levinson 2006). The corresponding event rate in a km^ detector is 
iV^ ~ lO^p yr~^ for a powerful blazar at a redshift z = I, and Nn ^ O.OS^p yr~^ 
for Mrk 421 and Mrk 501. In the case of microquasars up to a few events can 
be detected during a strong outburst if the viewing angle is sufficiently small 
(Levinson &: Waxman 2001; Distefano et al. 2002). The estimated neutrino flux 
from GRBs implies that only nearby sources can be individually detected by the 
upcoming experiments. However, the cumulative flux produced by the entire 
GRB population should be detectable assuming that cosmological GEBs are the 
sources of the observed UHECRs (Waxman Sz Bahcall, 1997). This requires that 
the rate of energy production of UHECRs in the energy interval 10^^ - lO^^ cV 
is comparable to the rate of 7-ray production by GRBs in the BATSE band, 
E ~ 10^^-^ erg Mpc-3 yr'^ (Waxman 1995). 
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